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ABSTRACT 

Double white dwarfs are expected to be a source of confusion-limited noise for the future gravita- 
tional wave observatory LISA. In a specific frequency range, this 'foreground noise' is predicted to 
rise above the instrumental noise and hinder the detection of other types of signals, e.g., gravitational 
waves arising from stellar mass objects inspiraling into massive black holes. In many previous studies 
only detached populations of compact object binaries have been considered in estimating the LISA 
gravitational wave foreground signal. Here, we investigate the infiuence of compact object detached 
and Roche-Lobe Overflow Galactic binaries on the shape and strength of the LISA signal. Since 
> 99% of remnant binaries which have orbital periods within the LISA sensitivity range are white 
dwarf binaries, we consider only these binaries when calculating the LISA signal. We find that the 
contribution of RLOF binaries to the foreground noise is negligible at low frequencies, but becomes 
significant at higher frequencies, pushing the frequency at which the foreground noise drops below 
the instrumental noise to > 6 mHz. We find that it is important to consider the population of mass 
transferring binaries in order to obtain an accurate assessment of the foreground noise on the LISA 
data stream. However, we estimate that there still exists a sizeable number (~ 11300) of Galactic 
double white dwarf binaries which will have a signal-to- noise ratio > 5, and thus will be potentially 
resolvable with LISA. We present the LISA gravitational wave signal from the Galactic population of 
white dwarf binaries, show the most important formation channels contributing to the LISA disc and 
bulge populations and discuss the implications of these new findings. 

Subject headings: binaries: close — stars: evolution, white dwarfs — gravitational waves — Galaxy: 
stellar content 



1. INTRODUCTION 

To obtain important information about a va- 
riety of astrophysical sources. Astrophysics must 
move toward observations beyond the electromag- 
netic spectrum. Several ground based gravitational 
radiation (GR) observatories are already working 
at full efficiency and collecting data (e.g. , LIGO, 
IThe LIGO Scientific Collaboration: B. Abbott! §007)). 
These detectors are sensitive in a high GR frequency 
regime (above a few tens of Hz), where signals are 
expected to come from mergers of compact stellar 
mass objects like neutron stars (NSs) and black holes 
(BHs). However, these types of events are infrequent, 
their rate estimates are bur dened with heavy uncertain- 
ties (e.g., Bclczvnski ct al. 120021: IKalogera et"!!] 120041: 
lO'Shaughnessv et al.. ,2008. ') and it is not clear whether 
any merger detections ar e likely given the current level 
of instrument sensitivity (|LIGO Scientific Collaboration! 
[200&) . The future generation of GR detectors in- 
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eludes LISA, the Laser Interfe rometer Space Antenna 
(|Danzmannlll996l : iHughesI 120061 and references therein), 
a joint ESA and NASA mission. It is a space-based 
all-sky interferometer which will be sensitive to lower 
frequencies (^ lO^'' — 0.1 Hz) than are ground-based 
GR detectors. The different frequency regime will al- 
low for the observations of sources of a different nature. 
The most promising sources for LISA are mergers of su- 
permassive and /or intermediate-mass black holes (i.e., 
ISesana et a l."2005V extreme mass ratio inspirals (EM- 
RIs; e.g., iGah 2001) of stellar -mass objects into super- 
massive black holes , and stellar mas s compact remnant 
binar ies (see e.g., iHils et ahl 119901 : iBenacauista et all 
[200l . 

At first, it was believed that contact W UMa sys- 
tems (.Mironovsk ii 1965) would dominate the GR sig- 
nal in the low fr equency range, but later work (e.g. , 
lEvans et al.lll987t ILipunov et allll987l: IHils et al.![l99a) 
has shown that close double white dwarf (WD) bina- 
ries are more important at low frequencies, and are a 

cniarant.ppd sniirce of crravitatirtTial wayeS ^Or LISA. The 
utB.edu.,, s.lars6n<§usu.ediL.gabriieI.l.wilIianjs@gmair.cQ™ tttt^ 

the Galaxy will put constraints on the star formation 
history of the Milky Way (MW), as well as the scale 
heigh t and shape of the thick disc and Galactic bulge 
( Ben acauista fc Hollev-Bockelmaiin| [2b06) . 

The Galactic double WDs will contribute to the LISA 
signal as both unresolved and resolved sources. At 
low frequencies (< 3 mHz; pertaining to binaries with 
orbital periods larger than 11 minutes), predicted bi- 
nary numbers are so high that most of the GR sig- 

^ For a circular binary, /gr [Hz] = 2/ (orbital period) [s]. 
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nals from individual systems will be unresolved, and 
they will form a confusion-limited 'foreground' noise. 
Within that frequency range, only GR sources that are 
relatively close to the detector or are strong GR emit- 
ters will stand out above the foreground noise. At 
higher frequencies (> 3 mHz; orbital periods < 11 min- 
utes), the numbers of double white dwarfs are rela- 
tively small and thus they are expected to be resolved, 
offering an opportunity to study uncertain parame- 
ters associated with binary evolution and WD physics 
(e.g., possible progenitors of SNe la and R CrB stars). 
LISA will be mostly sensitive to Galactic stellar-mass 
binaries, however a small contributio n due to extra- 
galactic systems was also predicted (iHils et al. 119901 : 
iKosenko fc Postnovll99llFarmer fc Phinnevl2003h . The 
gravitational wave radiation from coalescing double WDs 
was investiga ted using smoothed particle hydrodynam- 
ics by iLoren-Aguilar et all ()2005D . It is expected that 
the signal arising from merging double WDs would be 
detectable prior to the coalescence, though would not 
contribute significantly to the LISA gravitational wave 
signal overall, as the Galactic double WD merger event 
rate is predicted to b e as low as ~ 1 per century 
(jNelemans et al.ll2001bO . 

Among the objects which are expected to be im- 
portant verification sources for LISA are the com- 
pact AM Canum Ven aticorum (AM CVn) binaries (see 
IRoelofs et al.l (|2007al) for gravitational wave strain am- 
plitude estimates of five known AM CVn binaries). AM 
CVn systems are a sub-class of cataclysmic variables 
with < 1 hour orbits in which a WD accretes matter 
via Roche-Lobe overflow (RLOF) from what is believed 
to be a heli um-rich whi t e dwarf or a (semi-degenerate) 
helium star CSmakl llQGTt IWarnerifTQQSbt INelemans et al.l 
[2001a). These close binaries are faint and thus not read- 
ily detected electromagnetically, but are expected to be 
observable in gravitational waves (25 AM CVn systems 
have been confirmed to datef^. 

Besides double WD binaries, other close compact bina- 
ries includin g NSs and BHs have been studied b y other 
groups (e.g.. lNelemans et al.ll200Tbl: lCooravll2004D . Some 
few tens of these systems with orbital periods within the 
LISA sensitivity range are known in the Galaxy. There 
are also ~ 8 ultracompact X-ray binaries, which consist 
of a WD transferrin g matter through RLOF to a NS on 
~ 40 minute orbits (jWang fc Chakrabartvl 120041) . Dou- 
ble NSs are also known, however only one such system is 
observed with a period (of 2 4 hrs) within the LISA sen- 
sitivity range (jBurgav et al.l 120031 ). Close binaries with 
BHs (like BH-BH, BH-NS, BH-WD) have yet to be ob- 
served, but they are predicted to populate the Galaxy in 
smaller numbers than systems with NSs and/or WDs. 

Previous work on Galactic compact binary populations 
in context of LI SA has been don e with t he exclusion of 
RLOF svstems (IHils et all 119901: IPostnov & Prokhoro^^ 
119981: IBenacauista et al.l 120041: iTimpano et al.i i2006i) or 
with the entire RLOF - | - detached populati ons combined 
jNelemans et al. 2001b: Edl und et al. r2005l. In addition, 
iHils fc Bender (2000) studied AM CVn systems sepa- 
rately using an analytical approach, finding no significant 
increase in the confusion-limited noise for space-based 

http:/ /www. astro. ru.nl/'nelemans/dokuwiki/doku.phpTid 
=verification_binaries:am_cvn_stars 



GR detectors w i th the addition of RLOF binaries, while 
INelemans et al.l ()2004D calculated the GR amplitude of 
resolved WD binaries using population synthesis meth- 
ods (both detached and AM CVn systems), finding that 
LISA would resolve a total of 22000 Galactic double WDs 
(- 11000 detached, 11000 AM CVn). However, a more 
recent study which takes into account the local space 
density of AM CVn stars predicts that o nly ^10^ AM 
CVn b inaries will be resolved with LISA (jRoelofs et al.l 
I2007bf ). This estimate is still ~ 2 orders of magnitude 
above the current number of AM CVn systems which 
have so far been detected. 

In modeling our stellar population, we use an updated 
population synthesis method with recent results on mass 
transfer/accretion in compact object binaries, that is sig- 
nificantly different and more complex than the methods 
used in previous calculations. In this study, we assess 
the importance of detached and RLOF compact binary 
populations in the Galactic disc and bulge on the overall 
LISA sensitivity. In § 2 we describe our binary modeling 
techniques and signal calculations, in § 3 we show the 
physical properties (masses, orbital periods) of our dou- 
ble WD population, and show the most important evo- 
lutionary channels which contribute to the overall LISA 
signal arising from the four Galactic sub-populations (de- 
tached disc, RLOF disc, detached bulge, RLOF bulge). 
In § 4 we show the LISA GR signals, estimate the like- 
lihood of source resolvability, and comment on the de- 
tectability and selection effects of the potentially resolv- 
able systems within our Galaxy. In § 5 we compare our 
results with those of previous studies, and in § 6 we give 
a summary of our findings. 

2. MODEL DESCRIPTION 
2.1. Population synthesis model 

In this study, the combined Galactic population of 
RLOF and detached compact remnant binaries with or- 
bital periods spanning the LISA band is considered. We 
use the StarTrack population synthesis code described 
in detail in iBelczvnski et al.l (|2C)08l ) to evolve primor- 
dial binaries within the Galaxy from the Zero- Age Main 
Sequence (ZAMS) and calculate their properties. The 
code has undergone a nu mber of important updates since 
IBelczvnski et al.l ()2002f) . including detailed mass trans- 
fer calculations. Recent results on mass accumulation in 
context of ultracompact X-ray binary formation are given 
in Belczynski & Taam (2004)^hile for double WD bi- 
naries in iBelczvnski et al.i ( 20051 ) . 

We evolve the field population of single and bi- 
nary stars (50 % binarity) with solar-like metallicity 
{Z — 0.02). Single star evolution is followed from 
the ZAMS employing m odified analy tic formulae and 
evolutionary tracks from iHurlev et al.l (pOOQ) . For the 
evolution of binary stars, the calculations also include 
full orbital evolution with tidal interactions, magnetic 
braking, supernova natal kicks and g:ravita tional radi- 
ation emission (see IBelczvnski et al.ri2008L for formu- 
lae). ZAMS masses (Mzams) span the mass range 
0.08 — 100 M0. Single stars and binary primaries (Ma) 
are drawn from a 3-component broken power law ini- 
tial mass function (jKroupa et al.|[T993h . and secondary 
masses (Mb) are obtained from a fiat mass ratio distri- 
bution q =secondary /primary (e.g.. iMazeh et al..,1992,) . 
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We note as well that binary interactions leading to mass 
transfer events may alter the course of evolution for stars 
of a given initial mass. Initial orbital separations span 
a wide range up to 10^ R© and are drawn from a dis- 
tribution which is flat in the logarithm, while initial ec- 
centricities are taken fr om a thermal-equilibrium eccen- 
tricit y distribution (see lAbti 119831 : iDuauennov fc MavoH 
I199H) . Progenitors of ah binaries are initially formed 
on eccentric orbits. Tidal interactions between binary 
components circularize orbits before the first RLOF in 
a given progenitor system occurs. However, for a small 
fraction (10%) of double WD progenitors that are found 
initially on very eccentric orbits, the first RLOF is en- 
countered when the orbit is still eccentric. In such a 
case we assume an instant circularization at periastron 
(onew = cioid (1 " Sold) and Cnow = 0), and follow with the 
RLOF calculation. At the time of double WD formation 
all of the orbits are circular. Eccentric WD binaries are 
expected to arise from dynam ical interactions in globular 
clusters (jBenacauista II200TI ). and though we do not con- 
sider them here could provide a unique o pportunity for 
learning about WD structure with LISA (jWillems et al.l 

Hobs). 

For the population of Galactic disk binaries, we use 
a continuous star formation rate for 10 Gyr, while for 
the bulge population we use a constant star formation 
rate for the first Gyr, and no star formation thereafter 
(entire MW age is 10 Gyr). This translates to a present 
Galactic star formation rate of 4 Mq yr~^ (disc mass 
= 4 X 10^° Mq with a constant star formation rate for 10 
Gyr — 4 M© yr^-'^; see § 2.2), which is within reasonable 
agreement of the current Galactic star formation rate 
estimate of 3.6 M© yr~^ (|Cox 2000). We note however 
that it has been suggested that the star formation rate 
of the MW has been (mostly) decreasing exponentially 
with ti me, only reaching ^ 3.6 M© yr~^ at the current 
epoch (jNelemans et al.lfeoOlbl 12004 see sect. 2.2); the 
integrated mass in formed stars being closer to 8 x 10^" 
M©. Since the star formation history of the MW disc 
and bulge is not precisely known, we choose to model 
the Galaxy with simple star formation histories (1 Gyr- 
long 'continuous burst' (bulge), and constant for 10 Gyr 
(disc)), which to first order is a reasonable representation 
of the global star formation history of the Galaxy. 

The most uncertain phase in close binary evolu- 
tion affecting the orbital separation for low- and 
interme diate- mass stars is th e common envelope (CE) 
phase (jPostnov fc Yungelsonl [^006). Close binaries 
are expected to go through at least one common 
envelope event, and there are currently a few pre- 
scrip t ions of CE evolution in the literature (Wcbbinl^ 
T98l iNelemans fc Toutl[2005t Ivan der Sluvs et al... 20061 : 



Beer et al.l l2007| ) . In our simulations the CE phase is 



treate d using energy balance as discussed in iWebbinfl 
()1984[) . in which the orbital energy of the binary is di- 
minished at the expense of the unbinding of the donor 
envelope. Th e post-CE sep aration is governed by the 
parameters A (|de Kooll 119901 ) . a function of the donor's 
structure, and the highly uncertain a, the efficiency with 
which the binary orbital energy is used to unbind the en- 
velope. We have chosen to use a x A = 1. The effects of 
using other CE efficiencies have not been fully explored 
here, though lower CE efficiencies result in closer post- 



CE binaries, and a higher number of stellar mergers. 

At the current age of the Galaxy (10 Gyr) we extract 
all binaries containing two compact remnants: WDs, NSs 
and BHs in the gravitational frequency range: 0.0001 — 1 
Hz, which encompasses the LISA sensitivity range (or- 
bital periods of 5.6 hr - 2s). Henceforth when we refer 
to ^LISA binaries', we are referring to binary systems 
in our study within this GR frequency band. As pre- 
viously stated, we consider only white dwarf binaries 
when calculating the LISA signal. There are 5 types 
of WDs considered in the evolution: helium (He WD), 
formed by stripping the envelope off a Hertzsprung gap 
or red giant low-mass star; carbon-oxygen (CO WD), 
formed from progenitors with masses ^ 0.8 — 6.3 M©; 
oxygen- neon (ONe WD), formed from progenitors with 
masses ~ 6.3-8.0 M©; hybrid (Hyb WD), having a CO- 
core and a He-envelope, formed via the stripping of the 
envelope from a helium burning star, and hydrogen (H 
WD), formed by stripping the envelope from a very low 
mass (< 0.8M©) main sequence (MS) star. In our LISA 
GR calculations, we only consider the first four types of 
WDs, as hydrogen white dwarfs are more representative 
of brown dwarf-like objects and thus we will from now 
on refer to them as brown dwarfs. For the current study, 
the brown dwarfs (formed through binary evolution) do 
not play a large role in contributing to the LISA GR sig- 
nal due to their low mass (see § 5). We note however 
that in older stellar populations, the fraction of bina- 
ries which contain H WDs becomes more significant (see 
iRuiter et al.l l2009l for the calculation of the LISA sig- 
nal from MW halo double white dwarfs, which includes 
'hydrogen white dwarfs'). 

In compact binaries, gravitational wave emission is a 
strong source of angular momentum loss. The average 
rate of orbital angular moment um loss due to GR for 
e = binaries is calculated from iPetersI (|1964l ) 



dJg, -32 G^/^ M2 M2 ^Mp + 



dt 



(1) 



where G is the gravitational constant and c is the speed 
of light. At every time step in our calculations the evolv- 
ing system is checked for RLOF. Upon reaching contact 
(RLOF) , we assume that mass-transfer is GR-driven and 
the donor (Afdon) mass transfer rate is calculated via: 



Mrinn = A/don -D 



-1 dJgi / dt 



(2) 



where R LOF stabilit y is determined by the parameter D 
(see i.e.. lKing fc Ko lb 1995) 

n_5 1 l-/a (l-/a)(l+g)/3.nt+/a 

6 2 3(1 + q) q 

Jorb is the orbital angular momentum of the binary, /a is 
the fraction of transferred mass accreted by the WD of 
mass Mace (/a = 1 for stable RLOF), q = (A//acc/Mdon), 
and /3,nt = Ml^J{Mdon + M^^cf. The radius mass ex- 
ponent for the donor Cdon i s obtained from stellar mod- 
els in each time step (see IBelczvnski et al.l 120081 ) and 
is ~ —0.35 for WD donors in a phase of stable mass 
transfer. Stable RLOF between two WDs leads to an 
increase in orbital period where dynamically unstable 
RLOF leads to a merger. Sec iBelczynski et al.. (,2008 ) for 
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a more detailed description of treatment of mass trans- 
fer/accretion phases in StarTrack. 

2.2. Calibration and Spatial Distribution of Sources 

In our simulations both detached and RLOF systems 
have been distributed about the disc and bulge sepa- 
rately, assuming no correlation between position and age 
of the system. For the current study, we have neglected 
the halo population of double WDs, but t hese systems 
have been investigated in another work (jRuiter et al.l 
|2009( ). The density distribution of the Galactic disc is 
taken to have the following form: 



p(R,z) 



AttR^zo 



,-R/Rop\-z\/zo 



(4) 



in cylindrical coordinates, and the density distribution of 
the bulge has the form: 



inr^ 



(5) 



in spherical coordinates (bulge scale length rg = 

Nd and Nb represent the total numbers of simu- 
lated double WDs in the disc and bulge, respectively. 
We have calibrated our res ults using ste l lar G alactic 
disc and bulge masses from iKlvpin et al.l (|2002| ) , with 
masses of 4 x 10^" and 2 x 10^° M©, respectively. 
We have chosen i?o = 2500 pc and zq = 200 pc 
for the disc scale length and scale height while tq = 
500 pc with a r adial cut-off of 3500 pc for the bulge 
(jNelemansI I2003f) . While we note that even a scale 
height of 500 pc would be a reasonable choice for disc 
WDs (Majcwski & Siegcl 2002), we have chosen to use 
200 pc so that our results are more readily compara- 
ble to those of previous studies (iNelemans et al.ll2001bl : 
IBenacauista fc Hollev-Bockelmannll2006l ). 

Our original LISA double WD binaries which were 
born directly from StarTrack (8.4 x IC^) were the re- 
sult of the evolution of 20 x 10^ ZAMS binaries. Since 
evolving the whole Galactic population of LISA binaries 
is time-intensive, an interpolation scheme was developed 
in order to scale our number of StarTrack binaries to 
match those of the Milky Way bulge and disc by stellar 
mass. Probability distribution functions (PDFs) were 
constructed for different evolutionary channels as a func- 
tion of the primary and secondary masses and GR fre- 
quency. For RLOF binaries, it was only necessary to 
construct PDFs as a function of the masses since the or- 
bital period, and hence GR frequency, was then uniquely 
determined by the masses. We chose to construct PDFs 
(instead of simply scaling the original binaries by a fac- 
tor of > 300) so that our systems would span a range 
of frequencies and masses pertinent to said formation 
channel. The PDFs were constructed using a kernel den- 
sity estimator t hat best replica ted the characteristics of 
the population (| Williams! l2008[ ) . The KDE package for 
Matlab was used to generate the specific PDFs. This 
technique was particularly successful for well-populated 
evolutionary channels, where we could be sure that the 
range of masses and frequencies was well covered. For 
some of the least-populated channels, the full range of bi- 
nary masses and frequencies was not covered sufficiently 



to yield a smooth extrapolation of the data. Thus, appar- 
ent clumps of binaries that are comparable in number to 
our scaling factors are attributable to an under-sampling 
of the frequency and mass range for these rare evolution- 
ary channels. 

2.3. Lisa Signal Calculations 

The gravitational wave signal from the LISA data 
stream will be a time delay interferometry (TDI) variable 
in which the phases of the laser signals at each vertex of 
LISA are combined with the phases of other signals at 
delayed times in order to reduce the laser phase noise 
and accommodate the var ying armlengths of the con- 
stellation of the spacecraft (jRubbo et al.|[200^ . At low 
frequencies, the signal at any given vertex can be very 
well approximated by the Michelson signal: 



h{t) = -hab {ei£\ 



(6) 



where hat is the wave metric and £i and £2 are the 
unit vectors pointing along the two arms that join at 
the vertex. The transfer frequency is the frequency at 
which point the gravitational radiation wavelength be- 
comes comparable to the armlength of LISA. For high 
(e.g., above the transfer frequency) GR frequencies, the 
period of the gravitational wave is then less than the time 
it takes for light to propagate between LISA spacecraft 
detectors, and the low frequency approximation breaks 
down. The transfer frequency is given by /, — c/{2ttL) « 
0.01 Hz, where L is the LISA arm-length taken to be 
L — 5 X 10^ m. At frequencies above or near a more 
accurate representation of the LISA signa l is given b y the 
rigid adiabatic a pproximation jR ubbo ct _al ll200^ . but 
detailed analysis (jVecchio fc Wickhamii200^ have shown 
the low-frequency approximation should be adequate for 
most tasks at frequencies below ^ 30 mHz. In this work, 
we have used the rigid adiabatic approximation to calcu- 
late the LISA signal for sources with frequencies above 
0.003 Hz. 

Space-borne gravitational wave observatories like LISA 
are constantly in motion, changing their relative speed 
and aspect with respect to astrophysical sources on the 
sky. The sensitivity of the observatory to different grav- 
itational wave polarizations as a function of sky position 
is reflected in the antenna beam patterns, which have 
the highest gain in the direction perpendicular to the 
plane of the interferometer. As the observatory moves in 
its orbit, the gravitational wave signal is modulated in 
amplitude, frequency, and phase. Frequency (Doppler) 
modulation arises from the relative motion of the obser- 
vatory and the source, amplitude modulation arises from 
the sweep of the anisotropic antenna pattern on the sky, 
and phase modulation results from the detector's chang- 
ing response to the gravitational wave polarization state. 
The Michelson signal can be calculated using Eq. (|6]) by 
either holding the source fixed and putting all th e mo- 
tion of the detector into £1 and £2 (|Rubbo et al.l [2004') . 
or by holding the detector fixed and putting the motion 
into the source (Cutler 1998). Both approaches yield the 
same Michelson signal in the low-frequency limit: 



h{t) = ^A(t) cos 



2ti f {t')dt' + ipp{t) + if D{t) 



(7) 
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where A{t) is the ampUtude modulation, / — 2/Porb 
(where Porb is the binary orbital period) is the gravita- 
tional wave frequency for systems with zero eccentricity, 
and ifo is the initial phase of the wave at i = 0. The po- 
larization phase, (pp{t), represents the phase modulation. 
The Doppler phase, ipoit), gives the frequency modu- 
lation and the amplitude modulation is described by: 
A{t) = ^{A+F+)^ + (AxF^)2. The sensitivity factors, 
and are complicated functions of the position an- 
gles and orientation of t he binary, as well as the position 
of LISA in its orbit rsee lCutledlfOM IR.ubbo et al. l lMil 
for specific forms of these functions) . The plus and cross 
polarization amplitudes are given by: 

A+=2^>^^/3(^/)2/3(l + cos2^) (8) 

Ax=-4^X5/3(^/)2/3cos*, (9) 

where d is the distance to the binary, i is the an- 
gle of inclination, and the chirp mass is = 
( M^Msf^^ / (M u + Ms )^'^^ We have used the approach 
of iRubbo et al.l ()2004[ ) to calculate the timestreams, 
which are then added to produce the total observatory 
data stream, which is then Fourier transformed to pro- 
duce the illustrated frequency domain representation of 
the double WD Galactic foreground. 

The LISA instrument noise is simulated by assuming 
the power spectral density of t he noise i s ma de up of 
position (or shot) noise given bv lCornishI (|20Q1 : 

S-^p = 4.8 X lO"'*^ Hz-\ (10) 
and an acceleration noise (converted to strain) given by: 

^na = 2.3 X 10-« (^^^)' Hz-1. (11) 

These separate components are combined according to: 

Sn = 45„p + 85„a (1 + C0S2 (///,)) , (12) 

where /, is the transfer frequency. We roll off the accel- 
eration below /min = 10~^ Hz, so that 5'na(/ < /mill) = 

5'na(/min)- In reality, the LISA noise will probably not 
follow this simple power law all the way down to our 
choice of /min, but will begin to rise at a higher frequency 
below 0.1 mHz. 

3. RESULTS 

Out of our total population of compact remnant bina- 
ries with gravitational wave frequencies within the LISA 
sensitivity range, we note that ^ 76% are double WDs, 

24% are binaries involving a white dwarf and a brown 
dwarf (WD-BD), while other types of binaries make up 
the rest (i.e., < 0.5 % of binaries are NS-WD binaries; 
less than 0.1 % are NS-NS). Specific predictions for dou- 
ble compact objects wi th NSs and BHs are discussed in 
iBelczvns ki et all ([2008). 

Relevant for the current study of double white dwarfs, 
we find N = 34.2 x lO*" LISA binaries in the disc+bulge 
population (He, CO, ONe and hybrid WDsJE . Out of 

When binaries involving brown dwarfs are considered, A'^ = 
45 X 10^. 



this population iVrf = 24.8x10^ (73 %) and A^b = 9.4x10^ 
(27 %). We find 8.6 x 10^ (25 % of total) of disc binaries 
are detached and 16.2 x 10® (47 % of total) of disc binaries 
are RLOF, while for the bulge 3.9 x 10® (12 % of total) 
are detached and 5.4 x 10® (16 % of total) are RLOF. 
The total number of detached binaries among the total 
disc-Hbulge LISA population is 12.5 x 10® (37 %; 8.6 x 10® 
and 3.9 x 10® for the disc and bulge respectively) while for 
RLOF, often neglected in previous studies, the number 
is 21.6 X 10® (63 %; 16.2 x 10® and 5.4 x 10® in the disc 
and bulge, respectively). Within our RLOF population, 
> 99 % (21.4 X 10®) are AM CVn-like systems (a WD 
accreting from a helium or hybrid WD), the majority of 
which are COWD-HeWD binaries. 

The number of AM CVn systems which we find to 
currently exist in the MW can b e compared to the es- 
timate of INelemans et all ()2001aD . However, We remind 
the reader that in our population synthesis, although we 
do form AM CVn stars through the helium star channel, 
we only include binaries consisting of two degenera te ob- 
jects in the final results (see iNelemans et al.ll2001al . for a 
description of AM CVn formation channels). Thus, we 
cannot perfor m a direct comparison w i th the AM CVn 
popula tion of INelemans et all ()2001aD . INelemans et al.l 
(2001a) find the current number of Galactic AM CVn 
(double WD channel only) to range between (0.2 — 49) 
xlO®, depending on the assumed effective tidal coupling 
(see their § 3.4), and they find corresponding space den- 
sities of (0.4 - 1.7) xlO~'^ pc~^, respectively. Thus, we 
find that our AM CVn number estim ate (21.4 x 10®) 
is in re asonable agreement with that of INelemans et al] 
(j2001al) if in fact the double WD evolutionary pathway 
is an efficient channel for creating these systems. 

In the study of Roclofs ct al. (2007b), it was deter- 
mined that population synthesis studies were overesti- 
mating the space densities of AM CVn binaries, and 
that the true local space density is po = 1 ~ 3 x 10~® 
pc~'^, based on spectroscopic observations of 6 SDSS-I 
AM CVn stars. Adopting the disc-like density distribu- 
tion as discussed in § 2, we find that the local space 
density of AM CVn binaries from our population syn- 
thesis model is 2.3 X 10-5 pc^^fT^ which is stiU an order 
of magnitude above the most up-to-date observational 
estimate. It is becoming more clear that mass transfer 
between two WDs may not be as efficiently sustained as 
has been assumed in most population synthesis studies, 
and that a large fraction of these systems should merge 
upon r eaching contact ra ther than enter the AM CVn 
phase (jMarsh et al.ll2004D . 

3.1. Formation Channels 

We note that in our calculations, the total population 
of WD binaries with GR frequencies between 1 x 10^^ — 1 
Hz comprises roughly 6% of the total StarTrack Galactic 
population of double WDs. While we include all RLOF 
double WD systems that currently exist in the MW in 
our study, only 2.3% of Galactic detached double WD 
systems are accounted for, given the orbital period cut-off 
of ^ 5.6 hrs. In Tables 1, 2, 3 and 4, we show the evolu- 
tionary history (formation channels) of the most common 
LISA double WDs. We indicate the contribution (per- 

For a spherical volume with a radius of 200 pc. 
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centages) of said channel to the population of LISA WD 
binaries from that particular Galaxy component (disc, 
Tables 1 & 2, or bulge, Tables 3 & 4). The characters 
and numbers in parentheses in the central column of the 
tables represent the formation channel histories, and re- 
count various stages of stellar evolution of the primary 
and secondary star progenitors (see Table 1 caption for 
description) . 

In our formation channel notation, for a 'CO-He' bi- 
nary the CO white dwarf is the first formed WD, which is 
not necessarily the more massive WD (thus is not neces- 
sarily the primary star). In general, the detached systems 
have evolved from progenitors in which the stars initially 
have comparable massesO thus first RLOF mass trans- 
fer is more often dynamically stable and the binary does 
not undergo a CE at this stage (such as He-He detached 
double WDs). Typically the detached double WD pro- 
genitors go through only one CE event, where as in many 
cases the RLOF binaries undergo two CE events. 

We wish to point out that in the case of several RLOF 
binaries (e.g., those involving hybrid donors), there is 
a short-lived detached phase preceding the long-lived 
RLOF phase. For example, in the evolutionary history of 
a COHyb-Rl binary, the progenitor system only spends 
^ 10^ Myr as a detached CO-Hyb WD after the second 
CE phase before driven to contact via GR. Once contact 
is reached however, the binary will spend a long time 
(Gyr) as a stable RLOF system. This explains why the 
COHyb-Rl channel shows up in Table 2, but its pre- 
decessor COHyb detached channel is absent in Table 1 
(COHyb-Dl binaries only account for ~ 1% of the disc 
formation channels for LISA double WDs). 

In Figure [l] Figure [H Figure [3] and Figure |4] we show 
in (arbitrary) colour the number densities of the most 
common formation channels for the detached and RLOF 
populations in the disc and bulge. We do not include 
the less-populated formation channels in the plot or it 
would be very difficult to distinguish between channels. 
We do however include some formation channels involv- 
ing brown dwarfs, since even though these binaries do not 
significantly affect the overall GR signal (and we have ig- 
nored them in the signal calculation) , they are relatively 
abundant, specifically in the bulge population. 

It is expected that the disc would host binaries which 
have evolved from a wider variety of formation channels 
than the bulge, since the age of the disc binaries extends 
over a large range: ~ few hundred Myr to 9 Gyr. Note 
that high frequency systems only exist in the Galactic 
disc. There is a depletion of very high frequency (few 
minutes orbital period) RLOF systems in the bulge since 
the present bulge population only contains double white 
dwarfs with long-lived (> 9 Gyr) formation histories. All 
of the heavier (ONe-t-CO, CO-hCO, CO+hybrid) binaries 
originating from more massive progenitors which were 
born during star formation in the bulge have either long 
since merged, or in the case of RLOF binaries, the stars 
has long since reached contact and are now exchanging 
mass in RLOF on slowly expanding orbits0 Once any 
system encounters stable RLOF the binary's orbital pe- 

Many RLOF systems originate from systems in which the ini- 
tial primary star was much more massive than the initial secondary, 
leading to a CE. 

We also note that in certain cases, particularly for bulge 
RLOF systems such as some AM CVn binaries, the donor has 



riod increases as a consequence of conservation of an- 
gular momentum under continued mass exchange, since 
the less-massive (larger) WD is losing matter to the more 
massive WD. For the Galactic disc, some WD binaries 
are still relatively young, and are in an earlier phase of 
stable RLOF and thus we 'catch' these mass transferring 
systems at shorter orbital periods (GR frequencies above 
- 4.5 mHz; Porb < 7 - 8 min). 

Typical detached evolution: Detached double CO- WDs 
are formed through a variety of evolutionary channels. 
For the most prominent CO-CO channel (COCO-Dl, 
see Table 1), one particular example of evolution pro- 
ceeds as follows: Two MS stars (2.88 and 2.45 Mq) 
start out with a Porb of 14.2 days and an eccentricity 
of 0.54 when the MW is 9180 Myr old. At 9602 Myr, 
the more massive star (si) begins to evolve off of the 
MS and becomes a Hertzsprung Gap star. Shortly there- 
after si fills its Roche Lobe and begins to transfer mat- 
ter to the companion (interaction between component 
stars aids in circularizing the orbit; Porb ^ 4.5 days). 
RLOF is dynamically stable, and since si quickly be- 
comes the less-massive star during mass transfer Porb 
increases (to ^ 143 days), and RLOF stops at 9607 Myr 
when si is a red giant (0.44 M0; core mass 0.43 M©) and 
s2 (3.67 Mq) is still on the MS. si continues to evolve 
and shortly thereafter becomes a (naked) helium star 
(0.44 M0), as it has lost the remaining part of its de- 
pleted (by prior RLOF) envelope. At 9731 Myr, s2 has 
evolved off of the MS to become a Hertzsprung Gap star. 
By 9777 Myr, s2 has now evolved into an early AGB star 
(3.62 Mq), and the orbital period has decreased to 81 
days due to tidal interactions (expanding s2; si is still 
a helium star). s2 then fills its Roche Lobe, and since 
now the mass ratio is rather large (Mdon/-Wacc ~ 8) the 
mass transfer is dynamically unstable, leading to a CE 
phase. After the CE event the donor (s2) has become 
a Hertzsprung Gap helium star (0.82 Mq) and Porb has 
decreased by more than 2 orders of magnitude to 3 hours. 
The orbit then starts to decay slightly due to angular mo- 
mentum losses associated with GR emission, and there 
is a brief RLOF phase in which the Hertzsprung Gap he- 
lium star (s2) loses ~ 0.1 M©, half of which is gained 
by si (still a helium star, now 0.49 M©). When the 
MW is 9779 Myr old, s2 finally becomes a CO WD (pri- 
mary WD), si continues to burn helium and becomes a 
Hertzsprung Gap helium star at 9838 Myr (Porb = 2.7 
hours). Within 20 Myr, si then becomes a low- mass 
CO WD (secondary WD) finishing its helium star evolu- 
tion. At 10 Gyr (present) the system is found as a de- 
tached double CO WD with component masses of 0.70 
& 0.49 Mq, and Porb = 2.28 hours (/ = 0.24 mHz). The 
two WDs will be brought together in ^ 300 Myr due 
to GR emission. The system is expected to merge due 
to dynamically unstable mass transfer with a combined 
mass of 1.19 M© (pre-merger orbital period of ~ 1.3 min- 
utes). The total mass of the merged system does not 
exceed the canonical Chandrasekhar mass and so it is an 
unlikely Type la Supernova progenitor, though may re- 
sult in the form ation of a massive rapidly r otating WD 
(|Saio fc Nomot o 2004; Yoon fc Langej[2005[ ). 

Typical RLOF evolution: In Figure [5] we show a repre- 

reached a lower mass limit of 0.01 Mq at which point we no longer 
follow the evolution. 
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sentative evolution through RLOF for one of the Galactic 
double WDs, a COHe-Rl channel system (see Table 2). 
The primordial binary consists of two low-mass MS stars 
(2.60 and 1.40 M©) on a wide (Porb 31.0 year) and 
eccentric (e = 0.9) orbit, born when the MW is 5410 
Myr old. At 5966 Myr the more massive star (si) be- 
gins to evolve off of the MS. The orbit begins to shrink 
due to tidal interactions and by 6126 Myr, si is an early 
AGB star, the orbit has circularized (e — 0) and the pe- 
riod has decayed to 2.5 years. Within 1 Myr si becomes 
a late AGB star. Magnetic braking (convective enve- 
lope of si) dominates angular momentum loss and the 
orbit decays slightly, si fills its Roche-Lobe (6127 Myr) 
and mass transfer is unstable, leading to a CE phase 
and a drastic decrease in orbital period (from 2.4 years 
to 8.2 days), si then becomes a COWD (primary WD, 
0.64 Mq). At 8782 Myr, s2 evolves off of the MS. The or- 
bit decays due to magnetic braking (convective s2), and 
at 9085 Myr s2 initiates a CE phase resulting in another 
prominent period decrease (6.1 days to 30 mins) and the 
loss of its hydrogen-rich envelope. s2 then becomes a He 
WD (0.22 Mq). GR causes the WDs to come into con- 
tact a few Myr later (Porb ~ 3 mins), and the secondary 
(s2) fills its Roche-Lobe. The mass ratio is only ~ 0.34 
thus RLOF is stable, and an AM CVn system is born. 
Initially mass transfer is quite rapid (see Figure [5] top 
panel), but a slower mass transfer phase follows. The 
binary is found at 10 Gyr with Porb ~ 1 hour (^ 0.6 
mHz) and component masses of 0.83 and ^ 0.01 Mq 
for the CO and He WDs, respectively. Many of the 25 
observed AM CVn binaries' parameters are uncertain, 
and detection of these systems is biased against the typi- 
cal, longer-period systems with lower mass transfer rates. 
However this particular simulated binary's properties are 
somewhat similar to the very few long-period AM CVn 
binaries which have been detected (e.g., SDSS J1552 & 
CE 315), though the nature of the donors is still unknown 
(|Anderson et al.|[2005HRoelofs et al.ll2007bD ). 

3.2. Detached population characteristics 

Characteristic properties of the entire detached popu- 
lation are shown in Figure [6l Disc systems are found at 
all orbital periods (top panel) within the LISA sensitiv- 
ity range, with a drop in number at shorter periods due 
to the fact that once systems reach contact, double WDs 
either fall into the RLOF population or otherwise merge 
and drop out of the LISA population altogether. Sec- 
ondary (Ms, less massive WDs, middle panel) lie predom- 
inantly in four regions: '-^ 0.1 Mq (He WDs with CO pri- 
maries), a prominent peak near 0.3 Mq (mostly He with 
some hybrid WDs with mostly CO or He companions), a 
small clump near 0.5 Mq (CO secondaries with CO pri- 
maries) and a short peak at ^ 0.7 Mq (CO or sometimes 
ONe secondaries with either CO or ONe primaries). The 
He secondaries with very low mass (0.1 Mq) derive from 
the same formation channels as catacly smic variables 
(simil ar systems have been discussed in iPodsiadlowskil 
(|20O8')), in which a CO WD accretes from a low- mass 
main sequence star for several Gyr (note the higher num- 
ber of these systems in the Galactic bulge). Eventually 
the donor becomes fairly exhausted of its hydrogen sup- 
ply, and meanwhile has built up a significant helium core. 
The donor eventually is depleted of so much mass that 
it reaches the hydrogen burning limit (0.08 Mq), and 



thus becomes degenerate (in this case, helium rich; a he- 
lium WD). Many of the systems whose secondary mass 
peaks near 0.3 Mq are progenitors of AM CVn stars 
(COHe), while most of the binaries involving more mas- 
sive progenitors (CO, ONe) will likely merge once they 
reach contact. Primary WDs (Mp, more massive WD; 
bottom panel) display a peak at 0.3 — 0.4 Mq, though 
most primaries have masses 0.5 — 0.8 Mq, and a smaller 
number of systems are quite heavy (0.8 — 1.4 Mq). The 
low mass peak is comprised of He WD primaries, with He 
secondaries. The majority of the rest are CO WDs, with 
either He or CO secondary companions. In both the mid- 
dle and bottom panels, we note that there are relatively 
fewer heavy WDs in the bulge population. This is be- 
cause all of the more massive progenitors (e.g., CO-CO) 
have had a sufficiently long enough time to evolve off of 
the main sequence, encounter 1 or 2 CE phases, reach 
contact (GR is also more efficient at bringing together 
two CO WDs than two He WDs) and have either merged 
out of the population altogether or in some cases, have 
moved into the RLOF population. A two-dimensional 
greyshade plot is presented in Figure [71 where we show 
the relation between secondary and primary masses of 
LISA binaries in the Galactic disc, presented in terms of 
relative total percentages of LISA double WDs. 

3.3. RLOF population characteristics 

Characteristic properties of the entire RLOF popula- 
tion are shown in Figured The RLOF systems with the 
maximum Porb (lowest GR frequencies) appear at 75 
minutes (log(/gr)= —3.34). The shape of the period dis- 
tribution (top panel) is due to the fact that most sys- 
tems spend the majority of their time at periods of ^ 1 
hr during the dynamically stable RLOF phase. As indi- 
cated previously, 99 % of RLOF systems are AM CVn (a 
WD accreting from a He or hybrid WD). For most WD 
binaries, at the onset of stable RLOF orbital periods are 
on the order of only a few minutes, and mass transfer 
rates are initially high (see Figure [SI top panel). The 
donor quickly becomes exhausted of most of its mass, 
the mass transfer rate decreases and the period contin- 
ues to increase, though at a much slower rate. In the 
middle panel we show WD secondary masses, which is 
strongly peaked at ^ 0.01 Mq, with a slight drop-off to 
0.045 Mq. These WDs are comprised of mostly He WDs, 
some hybrid WDs and a small number of CO WDs. The 
'pile-up' of WDs at 0.01 Mq stems from the fact that 
once a AM CVn WD- WD progenitor reaches contact, 
initially the mass transfer rates are high (~ 10~^ Mq 
yr^^) and the donor is rapidly exhausted of mass, thus 
it is rare to catch AM CVn secondaries which are rel- 
atively massive (> 0.02 Mq). In the bottom panel we 
show the distribution of masses among the primary WDs, 
which are mostly CO WDs, with a small contribution of 
ONe WDs (masses > 1.3 Mq). We show in Figure [Hi a 
two-dimensional greyshade plot of RLOF LISA binaries 
in the Galactic disc, presented in terms of relative total 
percentages of LISA double WDs. 

4. LISA SIGNAL 

During a one-year observation period, the width of a 
resolvable frequency bin is A/ = 1 yr^^ w 3 x 10~^ Hz. 
Below 1 mHz, there are hundreds to thousands of binaries 
per resolvable frequency bin, causing the signal to be 
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confusion-limited. In our simulation, below ^ 0.45 mHz 
the only contributing binaries to the LISA signal are the 
detached double white dwarfs. 

The LISA spectra (amplitude densities, or spectral 
amplitudes hf) for the Galactic double WDs were 
calculated using the sop histicated LISA simulator of 
iBenacauista et all ()2004[) . and are shown in Figure [TUl 
Close binaries which evolve to contact, such as precur- 
sors to the AM CVn population or progenitors of white 
dwarf mergers, will spend a very short time (< 1 Myr) at 
GR frequencies > 0.01 Hz (Porb 3 minutes, provided 
they reach these short periods), and this is reflected in 
our original Star Track calculations. Thus, we artificially 
truncate the signal above 0.01 Hz since we simply do not 
have a large enough number of Star Track binaries above 
this frequency in order to accurately extrapolate them to 
the total number of Galactic systems here. However, we 
expect the number of LISA systems above 0.01 Hz to be 
< 10'^. In the bottom panel, we show the spectra from 
the four Galactic sub-populations, while in the top panel 
we show the full spectra for the combined population 
(grey) with a running median over 1000 frequency bins 
(white). The white curve represents our Galactic dou- 
ble WD foreground (see jj 4.2). We also show the forc- 
gro und estimates of lNelem ans et al.l (|2004D (dotted pink) 
and iHils fc Bended (2000) (dashed orange) for compari- 
son. The foreground of iNelemans et all ( 2004D was arti- 
ficially truncated beyond ~ 2 mHz, above which individ- 
ual binaries were expected to be resolvable in that work. 
It immediately comes to attention that the bulge bina- 
ries (green and mauve; lower panel) do not contribute to 
the GR signal above log(/gr) = —2.35 (/gr = 4.5 mHz) 
for reasons discussed in § 3.1. Also, the RLOF binaries 
do not contribute to the signal at low frequencies (be- 
low (log(/gi.)= —3.34, /gr = 0.45 mHz) since these are 
binaries which have been in a state of RLOF for a pro- 
longed period of time, and are on very slowly expanding 
orbits (largest orbital periods ~ 75 minutes). Addition- 
ally, once the donor mass drops below 0.01 Mq, we stop 
the calculation thus these very low-mass, low frequency 
systems do not contribute to the calculated LISA signal. 

The RLOF disc signal (blue; lower panel) becomes 
comparable to the detached disc signal (red; lower panel) 
at frequencies above ~ log(/gr)= —2.5 (/gr ~ 3 mHz). 
In our calculations, the AM CVn binaries spend only a 
short period of time (~ 10^ — 10^ Myr) as detached dou- 
ble white dwarfs, where as they may spend ~ Gyr in the 
RLOF phase. This is likely a consequence of the com- 
mon envelope prescription that we have employed, and 
most AM CVn sy stems go through a double CE (COHe- 
Rl channel). The iWebbinH ()1984[ ) CE prescription with 
a X A = 1 produces relatively close post-CE white dwarfs, 
and so it does not take long for the stars to be brought 
into contact via gravitational radiation once the detached 
double white dwarf has been formed. 

In the upper panel, we show our combined signal (grey) 
and the median foreground (white fine) alongside t he as - 
trop hysical foregrou nd estimates of lHils &: Bended (|2000f ) 
and INelemans et al. (2004) , both of which include de- 
tached and AM CVn lainaries. Note that our foreground 
is sig nificantly lower than the classic IHils fc Bendeii 
(|2000l) estimate at low frequencies, and is more compara- 
ble to the curve oflNelema ns et all (|2004[ ) in this regime 
(see below). Our foreground curve has been smoothed 



over 1000 bins, but some 'bumpiness' still exists at higher 
frequencies due to low number statistics. In contrast to 
some previous studies where the signal has been trun- 
cated at higher frequencies, we show the curve over a 
significant portion of the LISA bandwidth; up to 0.01 
Hz. Since our original StarTrack binary population was 
scarcely populated at high frequencies, consequentially 
we could not generate smooth formation channel PDFs 
(§ 2.2) for high-frequency systems. 

Our level of for eground signal is below that of 
INelemans etaU (|2004) for GR frequencies < 0.0013 Hz. 
This is in part due to the fact that at low frequencies, 
there is an overall lack of LISA binaries, since a large frac- 
tion of our LISA binary population are RLOF systems 
and hence do not exist at the low end of the frequency 
spectruml3 The median chirp mass of our double white 
dwarfs at frequencies below ^ 0.0013 Hz is ~ 0.38 Mq 
(disc) and 0.31 Mq (bulge). 

Despite the large number of RLOF systems in part of 
the low frequency, confusion-limited region of the Galac- 
tic gravitational wave spectrum (Figure [2] and Figure |4]), 
the amplitude of the gravitational wave signal at low 
frequencies is dominated by the detached double WDs. 
Mass-transferring binaries at the lower frequencies have 
undergone a significant amount of mass transfer, which 
leaves them with lower chirp masses than the detached 
binaries. As discussed previously, RLOF systems do 
not occur with orbital periods greater than ^ 2.5 hours 
(/gr < 0.45 mHz). Consequently the gravitational wave 
amplitude of the RLOF systems (both from the disc and 
from the bulge) does not contribute very strongly to 
the overall foreground noise in the low- frequency regime. 
The GR amplitude from the bulge is smaller than that 
of the disc for both (corresponding) detached and RLOF 
components of disc, primarily due to the fact that the 
stellar mass of the bulge is ~ 1/3 the mass of the disc, 
and contains fewer LISA double WDs. 

4.1. Transition Frequency 

The 'transition frequency' is said to correspond to the 
GR frequency at which the gravitational wave spectrum 
transitions from confusion-limited to individually resolv- 
able sources. The exact conditions which are necessary 
in order for this to occur (one possibility would be when 
the average number of binaries per frequency bin drops 
below 1) are not completely clear. The addition of RLOF 
systems to the detached population will increase the fre- 
quency at which the average number of binaries per bin 
drops below 1, simply by adding more sources to the 
population. This is particularly so since RLOF bina- 
ries are able maintain short orbital periods (for longer 
then their detached counterparts) once contact has been 
established. Furthermore, as the transition frequency in- 
creases, the spread of the signal due to the motion of 
LISA also increases, thus the signals will continue to 
overlap at higher frequencies than might be inferred from 
a simple analysis of the number of binaries per resolvable 
frequency bin. For a simple detached system in a circu- 
lar orbit whose orbital period is evolving solely due to 
the emission of GR, seven parameters are needed to fully 

^® We wish to point out that the number of binaries per resolv- 
able frequency bin at 0.0001 Hz in Nelemans et al. (2001) is a 
factor of ~ 2 higher than ours. 
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characterize the signaL Basic estimates frora information 
theory describe how the hmited amount of information 
in the LISA data stream can be mapped into an equiv- 
alent amount of information about detected binary sys- 
tems. A common rule of thumb is the "three-bin rule" : 
at least three (otherwise unpopulated) frequency bins 
of LISA data are needed to resolve a s even-parameter 
binary source (see iTimpano et alj 120061 and references 
therein for a discussion). One can reasonably expect that 
a mass-transferring system will require more parameters 
to characterize the GR signal, and consequently more 
frequency bins to individually resolve the source. There- 
fore, the average number of binaries per frequency bin 
is not a particularly good indicator of how the transi- 
tion frequency changes with the addition of the RLOF 
population. 

For the sake of comparison with previous studies, we 
have estimated the transition frequency using the 3-bin 
rule method. The point at which LISA binaries start to 
occupy one source per three frequency bins is found to 
be 1.5 mHz for t he detached LISA binar ies o nly, which 
is comparable to iNelemans et al.l (|2001bt ) and iNeleman^ 
(EoOBb) , who found a transitio n frequency of 1.5 — 2 mHz 
see also (jNelemans et al.ll2004l . for a discussion of the in- 
clusion of AM CVn binaries in the confusion limit cal- 
culation). However we find that when we consider the 
entire LISA population of double WDs, the transition 
frequency is increased to 2.4 mHz. We would like to 
however point out that these frequencies do not repre- 
sent the frequencies above which one expects to resolve 
all LISA binaries; there are still a number of populated 
bins at higher frequencies. 

The 3-bin rule method of estimating the frequency at 
which binaries start to become resolved is not a robust 
method, since it is expected that at these higher frequen- 
cies, the frequency evolution due to mass transferring bi- 
naries can contribute a significant additional spreading 
of the signal. Further, this method does not take into 
account important criteria such as Doppler, phase and 
amplitude modulation, and one should consider alterna- 
tive (more technical) methods in computing t he thr esh- 
old frequency region (e.g.. Chapter 2 of Ruite3[2009( ) . In 
the following subsection, we discuss an approach other 
than the 3-bin rule of thumb to identify the population 
of resolved sources. Instead of using a general estimator 
like the transition frequency, it is possible to use a pre- 
dicted signal-to-noise ratio for any particular system to 
dete rmine whether it is p otentially resolvable or not (see 
also ITimpano et"aIll2006D . 

4.2. Noise 

For LISA observations of length Tobs the sensitivity to 
a particular binary source can be expressed by construct- 
ing a simple estimator of the signal-to-noise ratio (SNR). 
For circularized compact binaries, long observations will 
narrow the bandwidth of the source, and the binary will 
appear in the Fourier spectrum as a narrow spectral fea- 
ture with root spectral den sity hb(f), given in terms of 
the strain amplitude ho as ()Larson et al.ll2000() 

hb{f) = ho^TZ . (13) 

LISA will detect binaries in the presence of compet- 
ing sources of noise. There will be two components to 



the noise in LISA analysis: instrumental noise in the de- 
tector itself, and irreducible astrophysical noise from the 
total population of Galactic close binaries. The instru- 
mental component of the n oise is given by the shape of 
the LISA sensitivity curve ([Larson et al ]f2000). In this 
work, the Galactic foreground noise is computed from 
the synthesized Galaxy itself (i.e., with StarTrack). The 
gravitational wave strength of every binary in the syn- 
thesized Galaxy is estimated from its parameters, and a 
running median of the total signal is computed as a func- 
tion of frequency. Since the total number of individually 
resolved binaries is expected to be small compared to the 
total population, the running median should be a good 
estimate of the residual astrophysical noise foreground 
that would result from a fully realized science analysis 
procedure. The foreground resulting from this procedure 
has been shown in Figure [10] (top panel). Additionally, 
we have calculated the median foreground signal which 
is expected to arise from the unresolvable sources alone 
(the 'Reduced Galaxy'), which is likely a more realistic 
estimate of the true Galactic foreground. In Figure [TTl 
we show the Galactic foreground alongside the Galac- 
tic foreground signal minus the signal arising from the 
sources which were determined to be resolvable using our 
signal-to-noise estimate technique, assuming a SNR ratio 
> 5, (see § 4.3). Alongside our Galactic foregrounds we 
show for comparison the standard LISA sensitivity curve 
as well as the average Michelson noise of LISA. The lat- 
ter curve (see § 2.3) is a more appropriate comparison of 
the LISA instrumental noise with our calculated gravi- 
tational wave signal from the Galactic binaries, as both 
curves (our signal and the Michelson noise) account for 
directional and freq uency dependencies of the wave as 
measured by LISA (iRubbo et al.[ l2004l ). The standard 
LISA sensitivity curve (Larson SJ) is often presented in 
the literature, and is an approximation of the spectral 
amplitude h f scaled from the barycentred sensitivity am- 
plitude ho from the Online Sensitivity Curve Generator. 

The ultimate level of the confusion foreground (and 
as a result, the ultimate SNR of any given source) will 
depend strongly on the actual data analysis algorithm 
used to identify and subtract sources. To date no fully 
realized implementation of such a procedure has been 
developed. The computed foreground shown here should 
be a good estimate to the output of a fully developed 
analysis procedure. 

4.3. Resolvable Sources 

Though millions of double WDs are expected to be de- 
tected with LISA, a much smaller number of these sys- 
tems are predicted to be resolved, depending upo n the 
physic al properties of each binary (see also Strocc r et al.l 
WM\ for a discussion). Due to the small s epara, - 
tions of these compact binaries, iCoorav et al.l ([20041 ) 
predicted that nearly 1/3 of the resolved LISA bina- 
ries will be eclipsing in the optical regime and thus will 
offer an opportunity for these systems to be studied 
with follow-up observations, allowing for the determi- 
nation of WD physical paramet ers (i.e., radii). It has 
been noted bv iNelemansI ([20091 ) however that the in- 
trinsic brightness o f white dwarfs was overestimated in 
ICoorav et al.l ([20041 ). thus the expected number of LISA 
double white dwarfs with electro-magnetic counterparts 
should be lower. .Nelemana (,2009. ) predict a smaller. 
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though still relatively optimistic number of LISA sys- 
tems which may have optical and/or near infrared coun- 
terparts; possibly as many as ^ 2000 sources (see their § 
3.3). 

The total effective noise in the detector (instrumental 
+ foreground) is given by the spectral amplitude hf{f), 
and the SNR is estimated as (|Larson "eFaLll2000l) 



SNR: 



hbif) ho\/Tobs 



hfif) hfif) 



(14) 



A binary is deemed "resolvable" (distinguishable from 
the confusion limited Galactic foreground) if its SNR is 
greater than some detection threshold (typical SNR > 
5). Additionally, resolved sources in the population are 
sorted into monochromatic sources and chirping sources. 
Monochromatic sources do not evolve appreciably over 
the LISA observing time, Tabs- A binary is considered 
to be chirping if its frequency changes by more than the 
LISA frequency resolution. A/ > 1/Tobs- This criteria 
is a conservative one, as the expectation is LISA will be 
able to detect frequency changes which are smaller than 
the frequency resolution A/ (Takahashi & Scto 2002). 
The detection of binaries with appreciable evolution in / 
(or 'chirp') is useful, since this frequency evolution allows 
for the determination of Ai, from which the distance to 
the so urce can be calculated (see e.g., iNelemans et al.l 

i2nnibD . 

The estimate of the SNR in Eq. (fT4|) will not be precise 
in the regime where orbital periods change on timescales 
short compared to the observation time, or where the 
shape of the instrumental noise changes appreciably as a 
function of frequency. However, this simple estimator is 
a good tool to use as a first cut in a search over a large 
population of sources where a calculation of the SNR for 
each binary would be computationally prohibitive. For 
the purposes of this paper, this estimator should be per- 
fectly adequate since in large part the binaries of interest 
are not chirping across a significant portion of the LISA 
band during Tabs- Once a set of Galactic binaries has 
been identified using the simple estimator in Eq. ()14|) . 
it is computationally plausible to consider a more ro- 
bust estimation of the SNR for each system which can 
then be constructed from the data analysis technique of 
choice. For example, with chirping binaries a good esti- 
mate of the SNR can be motivated from matched filtering 
techniques by integr ating over a binary's spectr al energy 
distribution, dE/df (|Flanagan fc Hughes! floM ). 

Using this simple SNR estimator to resolve LISA bi- 
naries, we find that a total of ~ 11300 double white 
dwarfs will be resolvable in our Galactic sample: ^ 600 
chirping and ~ 10700 monochromatic. This is inferred 
under an optimistic assumption that all binaries with 
a SNR > 5 have a chance to be resolved. Within the 
resolvable population, we find ^ 4000 AM CVn bina- 
ries a nd ^ 500 possib l e Dou b le Degene r ate S cenario - 
DDS; llben fc Tutukwl (Il98l : IWebbinkI (Il98l - SN la 
progenitors (detached COCO WD binaries which will 
merge within a Hubble time, whose total mass exceeds 
the Chandrasekhar mass limit). Out of these possible 
SN la progenitors, the majority are formed through the 
C0C0-D4 channel, which involves the most massive pri- 
maries on the ZAMS. As for the AM CVn resolvable bi- 
naries, 7 % are chirping and arise mostly from COHe-Rl 



and COHyb-Rl channels, in near-equal numbers. Both 
formation channels involve two CE events and thus if 
these systems will be resolved with LISA, they will serve 
as useful objects for understanding close binary evolu- 
tion. Additionally, a number of progenitor AM CVn bi- 
naries are resolvable (~ 3000; detached channels), but 
only ~ 100 of these are chirping. 

If we relax our criteria and assume that any binary 
with a SNR > 1 is potentially resolvable, we find ~ 
43000 resolvable double WDs (- 1000 chirping, - 42000 
monochromatic). Additionally, we find ~ 2800 LISA 
binaries with SNR > 10, ~ 300 of which are chirping 
sources. The binaries with large SNR (> 10) typically 
ah have orbital periods < 2000 s (/gr > 0.001 Hz) and 
have either just recently encountered RLOF, thus they 
still have relatively large chirp masses, or they are close 
to reaching contact. 

In Table 5 we show a breakdown of the 11300 (2800) 
potentially resolvable LISA double WDs which have SNR 
> 5 (SNR > 10), with percentages relative to the re- 
solvable population in question. The properties of re- 
solved systems reflect some of the selection effects that 
are expected to arise from using LISA to identify double 
WDs. Systems with short orbital periods (and conse- 
quently high gravitational wave frequencies) are favored 
both because the amplitude of the signal scales as /^^'^ 
and because the number of binaries per bin drops with 
increasing frequency. Furthermore, systems with large 
chirp masses are favored simply because the amplitude 
of their GR scales as A4^^^. Thus, we see that the re- 
solvable population is biased toward high mass, high fre- 
quency (short period) systems. Realistically speaking, it 
is still unclear as to which method(s) of detecting resolv- 
able binaries with LISA will be most successful, and so 
we have discussed the likelihood of potentially resolving 
individual binaries using a few different methods. 

Detached resolved sources. A large number of the re- 
solved sources arc of the HeHe-Dl channel, where two 
initially low- mass (~ 1 Mq) evolved stars lose their en- 
velopes during mass transfer phases. HeHe-Dl are not 
only the most numerous channel of the detached sys- 
tems (Table 1 and Table 3), but also extend to high 
GR frequencies. Lumping all of the COCO-D forma- 
tion channels together, COCO detached binaries make 
up ~ 13 — 18 % of the resolved population, and poten- 
tial DDS progenitors make up 5% of the total resolved 
population. Though HeCO-Dl double WDs are more 
numerous compared to COHe-Dl binaries, fewer of them 
are resolved. We note that COHe-Dl progenitors involve 
two CE phases (HeCO-Dl only one), and overall occupy 
closer orbits thus contributing to higher GR frequencies. 

RLOF resolved sources. For RLOF systems, the ma- 
jority are of AM CVn-type. The most common RLOF 
channel is COHe-Rl, followed by COCO-Rl. COCO- 
Rl binaries are not common, but they lie on very close 
orbits, consisting of a massive 1 M0) CO WD ac- 
creting matter from a low-mass (< 0.2 Mq) CO WD. 
The resolved COHe-Rl AM CVn systems account for 
most of the total number of resolved AM CVn bina- 
ries. Since the formation of these systems involves two 
CE events, future LISA observations may put some very 
useful constraints on this rather uncertain evolutionary 
phase which is crucial to the formation of close binaries. 
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However, when comparing relative formation scenarios of 
AM CVn binaries, we must keep in mind that the pre- 
cise behaviour of the stars upon reaching contact - and 
the exact physical conditions of the donor at contact - 
are stih unclear, as are the physical characteristics of 
the AM CVn progenitors themselves. Further observa- 
tions and detailed modeling will be an important step to- 
ward understa nding the donors in AM CVn type systems 
(jPelove et alj 12007V For example, deviations from the 
equilibrium evolution of AM CVn systems, e.g., brought 
on upon by changes in the mass transfer rate following a 
nova explosion, result in a measurable frequency evolu- 
tion which can greatly differ from that which is expected 
for the long-term, equilibrium evolutio n case. This was 
already shown by iStroeer fc NelemansI ()2009l ) ; if short- 
term changes in the evolution of AM CVn systems are 
neglected, it can lead to erroneous estimations of some 
of the binary parameters. 

5. COMPARISON WITH PREVIOUS STUDIES 

The number of binaries import ant in con t ext o f the 
LISA GR signal was estimated by iHils et aH ()1990( ) and 
in their work they predicted 3 x 10'' Galactic close dou- 
ble WDs (orbital periods < 1 day) based upon a sur- 
face density star formation rate. However that estimate 
seemed to over-predict the observed local density of dou- 
ble WDs available at that time. They concluded that 
even if they decreased the estimated space density by a 
factor of 10 (decreasing the total number of Galactic dou- 
ble WDs), close WD binaries would cont inue to dominate 
the LISA GR signal at low frequencies (IHils et al .11 19901 
^ 1). Population synthesis studies fe.g.. INelemans et al.l 
\20m which incorporate current estimates of Galactic 
star formation rates and SNe rates result in a higher 
number of Galactic double WDs (> 21 x 10^0, albeit 
the majority of these would have orbital periods too large 
to be detectable with LISA. Still, we note here that cur- 
rent estimates of the number of double WDs with GR 
frequencies detectable with LISA are a factor of ^ 10 
higher than those used in 1990. 

To populate the Galaxy (the disc only), IHils et ahl 
()1990[) used the density distribution 

p = poe-^/^°e"l"l/"°, (15) 

where po is the central density, R and z are galactocen- 
tric cylindrical coordinates, and Rq and zq are the ra- 
dial scale length a nd vertical s cale h eight of the Galactic 
disc, respectively. IHils et ahl ()1990D used zq = 90 pc in 
their study, and while though a reasonable assumption 
for the scale height of the thin disc, likely led to their 
(too) low estimate of double WDs. More recently, it 
has been noted (e.g.. iBenacquist a fc Hollev-Bo ckelmand 
120061 that the scale height for the binary popu lation may 
be mu ch higher (zq ^ 200 pc) than that used in lHils et al.l 
(fT990l) . 

The IHils k Bende r (2000) curve is a modified result 
based on IHils et al.l (11990), who in calculating the GR 
amplitude in 1990, used too low of a space-d ensity esti- 
mate for double WDs. In the original study of Hil s et al.l 
()1990[) which neglected RLOF systems. Galactic dou- 
ble WDs begin to dominate the GR amplitude at ~ 0.1 

We find ~ 5 X 10** double WDs in the Milky Way. 



mHz, and continue to dominate until ~ 1.6 mHz, where 
the average number of sources per resolvable frequency 
bin drops below 1 (similar results for the transition fre- 
quency were found in their st udy that included RLOF 
systems; 'Hils fc BendeH (pOOOf )). The analytical study of 
Hils et al., (il990() was not directly comparable to more 
recent studies of the LISA fore ground which u tilized 
more sophisticated techniques. iTimpano et al.l (I2006D 
perfor med a re-calculation of the analvtical iHils et al.l 
(|1990() GR foreground for Galactic binaries in or der to 
compar e the results directly to th ose of Nclcmans et al.l 
(|2001bD f see ITimpano et al.l [20061 their Figure 10). Lit- 
tle was known about the true spac e density of dou - 
ble white dwarfs at the time of the | Hils et al.l (|1990[ ) 
study, and since ITimpano etahl (|2006i r was based on 
these results which underestimated the space density 
(and hence total number) of Galactic white dwarf bi- 
naries, the ITimpano et"aLl (|200g ) study had a 90% de- 
crease in source densit y as compared to the study of 
INelemans et al.l (|2001bD . This reduction in number of 
systems resulted in a decrease in the foreground GR am- 
plitude. Until more recent observations of double white 
dwar fs were available to better constrai n the space den- 
sity (IWarnerill995t IRoelofs et al.ll2007bD . the error went 
undetected for some time, as this decrease was serendip- 
itously counteracted by an increase in amplitude im- 
posed b^;_thfi_J_too}_Jarge double WD chirp masses used 
in the IHils et all (|1990t) an dljim pano et al. (2006) stud- 
ies. In s hort, the s t udies of ITimp ano ct al. (2006) (recal- 
culated IHils et all (| 19901) 1 and [Nclcmans et al. ( 20019 ) 
have comparable foreground levels, albeit for different 
physical reasons. An increased number of observa- 
tions of WD binaries since the time of the H ils ct aL| 
C1990 !) study (i . e., SP Y project [NaBiwotzki et al.li200l 
Ande rson et al.l 120051) have confirmed that th e (lower) 
WD masses used by INelemans et al.l (|2001b[ ) coincide 
more closely with th e masses o f WDs in binaries, where 
the masses used in IHils et al.l (I1990D do not (see also 
iNelemansI [2003H ITimpano et al.l 120061 for further dis- 
cussion). 

Based on lEvans et"all ()1987D . INelemans et al.l (|2001bf ) 
used the criterion that the transition frequency occurs 
when the average number of binaries per bin drops below 
1, and they found the transition at ~ 1.6 mHz. They also 
found that 12124 detached double WDs will be resolved 
with LISA, assuming all systems above the LISA sensi- 
tivity limit (a signal-to-noise of 1) are det ectable. In their 
study which included AM CVn svstems. INelemans et al.l 
(|2004f ) found a transition frequency closer to ~ 2 mHz, 
and that a total of ~ 22000 Galactic double WDs (half 
detached and half RLOF) would be resolvable. For a 
signal-to- noise ratio of 1, we find that ~ 43000 double 
white dwarfs with GR frequencies < 0.01 Hz arc poten- 
tially resolvable with LISA. Ho wever, in order to perform 
a reasonable comparison with INelemans et all (|2001bD , 
we need to consider only those binaries whose GR fre- 
quencies are above ~ 2 mHz (and below 0.01 Hz; see 
§ 4). In this case, we find that ~ 25000 double white 
dwarfs may be resolvable with LISA. 

In previous studies, it was assumed that the level of 
the Galactic foreground will not contribute above the fre- 
quency at which the average number of sources per bin 
drops below 1 (coined the 'transition frequency'). Since 
those studies do not calculate the GR signal beyond the 
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transition frequency, above ~ 2 mHz their results cannot 
be compared with ours. We find that the frequency at 
which the WD foreground drops below the instrumental 
noise level (~ 6 — 8 mHz) is above the transition fre- 
quency which was found in previous studies. In contrast 
to previous work, it is expected that at higher frequen- 
cies (> 3 mHz) the RLOF systems may contribute more 
significantly to the LISA signal than detached systems. 
This has implications for the 'transition frequency' be- 
tween unresolved and resolved systems: the transition 
frequency is shifted to higher frequencies with the addi- 
tion of RLOF systems. 

Since we do not follow the evolution of non-degenerate 
stars in our LISA calculations, the binaries involving 
'pre'-hydrogen WD donors (i.e., CVs) are not shown on 
the figures, though we note they would make some contri- 
bution to the LISA data stream. We do not include the 
'brown dwarfs' in our LISA signal populations, though 
we note that binaries involving these brown dwarfs will 
not contribute significantly to the Galactic gravitational 
wave background as these binaries have very low chirp 
masses. Since the formation of such objects takes on 
average several Gyr, it is however important to include 
them when consider ing very old stellar populations such 
as the Galactic halo (jRuiter et al.ll2009l ). where these sys- 
tems comprise over 70 % of the entire remnant binary 
population. In that work, it was found that the LISA 
GR signal from halo white dwarfs will be at least an or- 
der of magnitude below the signal arising from Galactic 
(disc-|-bulge) white dwarf binaries. 

6. SUMMARY AND CONCLUSIONS 

We have used the population synthesis binary evolu- 
tion code Star Track and the detailed LISA simulator of 
[Benacquista et al. (2004) to calculate the gravitational 
wave foreground arising from close white dwarf binaries 
in the Milky Way disc and bulge. Our population syn- 
thesis includes both detached and mass-transferring bi- 
naries, and it is found that the mass-transferring (RLOF) 
binaries begin to contribute significantly to the fore- 
ground noise at high GR frequencies (> 3 mHz, Figure 
10 bottom panel). 

Our results are in general agreem ent with the analyt- 
ical study of iHils & : Bende rl (120001) and the population 
synthesis studv oflNclcmans et al.l ()2004D for the low fre- 
quency regime: the LISA signal is dominated by de- 
tached binaries. However, in contrast to previous stud- 
ies, it is found that, despite the fact that the 'transition 
frequency' is found to be ~ 2.4 mHz using the 3-bin 
rule of thumb, at frequencies above ^ 3 mHz the pres- 
ence of RLOF binaries will likely make it difficult to re- 
solve individual sources, and we find that the Galactic 
foreground dominates the LISA noise up to at least ~ 6 
mHz. This is because the RLOF binaries add to the to- 
tal number of LISA systems at high frequencies, and at 
high frequencies the RLOF systems chirp backward and 
spread out over more frequency bins, which offers a new 
set of challenges for LISA data analysts. The details of 
the true orbital period evoluti on due to mass transfer 
and possible tidal e ffects (e.g.. IDelove &: BildstenI 120031 : 
iBildsten et al.l [20061) are not well understood, and this 
uncertainty may impact the ability of future data anal- 
ysis schemes to successfully resolve these systems. It is 
starting to become clear that such effects should be taken 



into acco unt when studying the population of compact 
binaries ([Racine et al.l 120071 ). One new and interesting 
result of our study is the possibility that RLOF sys- 
tems which have recently entered the mass transfer phase 
(small periods, higher chirp masses) can add to the con- 
fusion foreground at higher frequencies than were previ- 
ously predicted for the double white dwarf population. 

It should be repeated here that in this study, we em- 
ployed one model prescription of CE evolution. This 
important phase of binary evolution may have a large 
impact on a synthesized stellar population, and various 
model prescriptions of CE may lead to a different number 
of Galactic binaries and binary physical properties, which 
in turn affect the shape and strength of the LISA GR sig- 
nal. The number of close {LISA) binaries could poten- 
tially be decreased if the CE efficiency is low (more merg- 
ers) or if angular momentum balance rather than energy 
ba lance is assumed du ring the CE phase. It was shown 
in iRuiter et al.l ([20061 ) that when a n angular momen- 
tum balance prescription is adopted ([Nelemans fc Tou3 
120051 ). post-CE orbital separations are overall larger, and 
fewer binaries will enter a CV-like RLOF phase between 
a white dwarf and a main sequence star, which is a cru- 
cial step for some of the older (bulge C0He-R2) LISA 
binaries. 

We have shown that there is a large population of po- 
tentially resolvable binaries in the Galaxy. The approach 
of estimating resolvable systems described in § 4.3 ('SNR 
> 5'), relies on assuming that signals arising from indi- 
vidual binaries are separable, i.e., can be demodu lated 
using some sort of matched filtering technique (see iSetol 
[20021 and references their-in). If detected, once their 
properties are measured these systems may aid to solve 
several remaining problems in binary evolution theory. 
Especially interesting would be the detection of any dou- 
ble white dwarf binary with a mass close to or over the 
Chandrasekhar mass (potential SN la progenitor or pos- 
sible pre-NS accretion induced coUapse object), or the 
measurement of the orbital periods of systems which have 
just emerged from a common envelope (the phase that 
determines the formation of most compact object bina- 
ries). However, we reiterate that the detection of resolv- 
able binaries will depend on the development of sensitive 
detection schemes and data analysis techniques. 
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TABLE 1 

Detached Disc double WD Formation Channels'^ 



Channel'' 


Evolutionary History*" 


rcl. to Disc*"^ 


HeHe-Dl 


MT1(3-1) CE2(10-3;10-10) 


8% 


HeCO-Dl 


MTl(3-l/2/3) CE2(10-5;10-8) 


6% 


COCO-Dl 


MTl(2/3/4-l) CE2(4/7-5;7-8) MT2(7-8) 


7% 


COCO-D2 


MTl(2/3/4-l) CE2(ll-5;ll-8) 


3% 


COCO-D3 


MTl(2/3/4-l/2/3) CE2(7-5;7-8) 


1% 


COCO-D4 


MTl(2/3-l) MTl(8/9-l) CE2(ll-3;ll-7) MT2(ll-8) 


1% 


COCO-D5 


MTl(3/4-l) CE2(7-3;7-7) MT2(7-8) MT1(8-11) 


< 1% 


COHe-Dl 


CE1(6-1;11-1) CE2(11-3;11-10) 


5% 


the rest 




4% 


total 


8.6 X 10'= 


35% 



GR freq. range 1 X lO"* - 0.01 Hz (5.6 hrs -200 s orbital periods). 
^ In channel names XXYY, XX represents the first-formed WD. 

CE - common envelope phase; MT - mass transfer (RLOF) phase. Numbers following CE or MT indicate the binary 
component donor star, either 1 (initially more massive star on the zero-age MS [ZAMS]) or 2 (initially less massive 
star). The numbers in parenthesis correspond to the donor and accretor types; The foUowring can be used as a guide 
HBelczvnski et al.|[200D : 0: MS star < 0.7 Mq 1: MS star > 0.7 Mq 2: Hertzsprung gap 3: red giant 4: core helium 
burning 5: early AGB 6: thermally pulsating AGB 7: MS naked helium star 8: Hertzsprung gap naked helium star 
9: giant branch naked helium star 10: helium WD 11: carbon-oxygen WD 12: oxygen-neon WD 16: hydrogen WD 
17: hybrid WD. In the CE phase, the stars prior to and after the semi-colon represent the stages of evolution at the onset 
of and following the CE event, respectively. 

Percentage is relative to LISA binaries in said component (Disc or Bulge). 



TABLE 2 

RLOF DISC DOUBLE WD Formation Channels'* 



Channel 


Evolutionary History 


rcl. to Disc 


COHe-Rl 


CEl(5/6-l;8/ll-l) CE2(11-3;11-10) MT2(11-10) 


47% 


COHyb-Rl 


CEl(5/6-l;8/ll-l) CE2(ll-3;ll-7/17) MT2(ll-7/17) 


8% 


HeCO-Rl 


MTl(2/3-l) CE2(10-5;10-8) MTl(lO-ll) 


5% 


the rest 




5% 


total 


16.2 X 10» 


65% 



Tablenotes are the same as in Table 1. 
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TABLE 3 

Detached bulge double WD Formation Channels*' 



Channel 


Evolutionary History 


rel. to Bulge 


HcHc-Dl 
COHc-Dl 
HeCO-Dl 
COCO-D 
the rest 


MT1(3-1) CE2(10-3;10-10) 
CE1(6-1:11-1) CE2(11-3;11-10) 
MTl(3-l/2/3) CE2(10-5;10-8) 
various channels 


20% 

9% 
5% 
3% 
5% 


total 


3.9 X 10" 


42% 




* Tablenotes are the same as 


in Table 1. 

TABLE 4 

RLOF Bulge double WD Formation Channels^ 




Channel 


Evolutionary History 


rel. to Bulge 


COHe-Rl 
COHe-R2 
the rest 


CEl(5/6-l;8/ll-l) CE2(11-3;11-10) MT2(11-10) 
CE1(5-1;8-1) MT1(8-1) MT2(ll-l/0) MT2(11-10) 


37% 
20% 
1% 


total 


5.4 X 10** 


58% 



* Tablenotes are the same as in Table 1. 
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TABLE 5 

Expected number of resolvable LISA double WDs 



Type'' % SNR > 10 % SNR > 5 



Detached 

HcHc-Dl 8.4 9.4 

HcCO-Dl 6.3 5.5 

COCO-Dl 1.2 1.8 

COCO-D2 3.6 3.2 

COCO-D3 2.3 1.8 

COCO-D4 9.1 5.3 

COCO-D5 1.4 1.4 

COHe-Dl 23.0 22.3 

COHe-D2 < 1 

COHyb-Dl 1.2 1.2 

the rest < 1 < 1 

total 57.0 53.0 
RLOF 

COHc-Rl 17.3 21.4 

COCO-Rl 10.2 11.0 

COHyb-Rl 9.9 8.4 

HcCO-Rl 2.6 2.6 

HcHc-Rl 1.6 1.9 

HybCO-Rl 1.2 1.5 

the rest < 1 < 1 

total 43.0 47.0 



TOTAL 2800 11300 



Formation channels producing resolvable LISA double WDs in the Galaxy, assuming that every binary with SNR > 10 
or SNR > 5 with a gravitational wave frequency below 0.01 Hz has a chance to be resolved. Relative percentages (with 
respect to the total resolved population for the appropriate SNR threshold) for each channel are indicated. 
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O.OOOI 0.001 0.01 

Fig. 1. — Number density (n = (dN/df); A / is the size of a resolvable bin for a 1 year observation time, 1/Tobs = 30 nHz) of the 
most prominent detached LISA double WD evolutionary channels for the Galactic disc. Some channels from Table 1 have been left out 
for clarity, and we additionally show some of the more prominent LISA binary channels involving brown dwarfs (neglected in the signal 
calculations). Note that as the frequency increases, so does the concentration of resolvable bins. 
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O.OOOI 0.001 0.01 

Fig. 2. — Number density of the most prominent RLOF LISA double WD evolutionary channels for the Galactic disc. All binaries with 
GR frequencies < 0.45 mHz involve binaries with brown dwarf donors. 
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Fig. 3. — Number density of the most prominent detached LISA double WD evolutionary channels for the Galactic bulge. 
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Fig. 4. — Number density of the most prominent RLOF LISA double WD evolutionary channels for the Galactic bulge. Note that 
COHe-Rl is descended from COHe-Dl in Figure |3] in which the binary undergoes two CE phases. 
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Fig. 5. — Characteristic evolution as described in § 3.1 of a RLOF system with a CO WD accretor and a He WD donor, from the AM 
CVn formation channel COHe-Rl starting at the ZAMS. Top panel shows the mass transfer rate evolution (with the critical Eddington 
rate shown), middle panel shows the mass change of the two components, while the bottom panel illustrates the orbital period evolution. 
The two CE phases coincide with the prominent drops in orbital period when the donor (initial primary) is a late AGB star, and then 
when the donor (initial secondary) is a red giant. Stable RLOF begins at the period minimum of ~ 3 min, when the Galaxy is 9.1 Gyr old. 
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Fig. 6. — Physical properties of the Galactic population of deta<;hed LISA double WDs. The solid line shows the disc binaries, where as 
the dotted line shows the contribution from the bulge. Bottom and middle panels show the distribution of primary and secondary mass 
(more massive and less massive WD, respectively), while the top panel shows the orbital period distribution. Note the different bin size 
and scales used on the axes. 
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Fig. 7. — Primary and secondary WD masses for Galactic disc detached LISA binaries (corresponding to the solid line histogram in 
Figure [6J. The greyscale shows the relative contributions (percentage wise) relative to the Galactic disc LISA double WD population. 
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Same as Figure |6] but for RLOF systems. 
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Fig. 9. — Primary and secondary WD masses for Galactic disc RLOF LISA binaries (corresponding to the solid line histogram in Fig 
The greyscale shows the relative contributions (percentage wise). Note the different scaling on the y-axis as compared to Figure [7] 
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Fig. 10. — Spectra (amplitude density) of the gravitational wave signal for LISA galactic WD binaries. Bottom Panel: Disc detached 
(red), bulge detached (mauve), disc RLOF (blue) and bulge RLOF (green). Top Panel: The combined population is shown with no 
smoothing (grey), and with a running median over 1000 bins (white). Additionally, the foreground estimates of Nclemans et al. (200j) 
(dotted pink) and Hils_ls JSeildeEl (.2000) (dashed orange) are shown for comparison. Both Panels: Note that the signals have been truncated 
at frequencies above 0.01 Hz (see text). The standard LISA sensitivity curve (black long-dash) for a signal to noise ratio of 1 is also shown. 
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Fig. 11. — Median spectra shown in Figure 10 (smoothed): the median Full Galaxy alongside the median Reduced Galaxy (red; full 
median signal with resolved sources removed), the standard LISA sensitivity curve (dashed line; SNR=1, also shown in Figure 10) and the 
simulated LISA Michelson noise curve (§ 2.3). As in Figure 10, the spectra have been truncated above 0.01 Hz. Both the Full Galaxy LISA 
double white dwarf foreground (blue line) and the Reduced Galaxy foreground curve (red) fall below the Michelson noise at roughly ~ 8 
mHz, where the curves are relatively noisy. When one considers the median foreground spectra against the standard sensitivity curve, the 
Full Galaxy foreground drops below the sensitivity curve at ~ 7.5 mHz, while the Reduced Galaxy foreground drops below the sensitivity 
curve at '~ 6 mHz. The (black) noise curves differ in sensitivity level and shape since the underlying noise levels used to generate each noise 
curve are different. In our study we have used the Michelson noise curve to represent the LISA instrumental noise; the standard sensitivity 
curve is shown for comparison. 



